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We use our current understanding of the all-order singularity structure of gauge theory amplitudes 
to probe their high-energy limit. Our starting point is the dipole formula, a compact ansatz for the 
soft anomalous dimension matrix of massless multi-particle amplitudes. In the high-energy limit, 
we find a simple and general expression for the infrared factor generating all soft and collinear 
singularities of the amplitude, which is valid to leading power in \t\/s and to all logarithmic 
orders. This leads to a direct and general proof of leading-logarithmic Reggeization for infrared 
divergent contributions to the amplitude. Furthermore, we can prove explicitly that the simplest 
form of Reggeization, based on the absence of Regge cuts in the complex angular momentum 
plane, breaks down at the NNLL level. Finally, we note that the known features of the high- 
energy limit can be used to constrain possible corrections to the dipole formula, starting at the 
three-loop order 
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1. The high-energy Umit of gauge ampUtudes 

The high-energy Umit of scattering ampUtudes, defined as the limit in which the center-of- 
mass energy of the colhding particles, -y/^, is much larger than other relevant kinematic invariants, 
has been an interesting and lively subject of studies for several decades [1, 2, 3]. In renormalizable 
field theories, the high-energy limit is characterized by the presence of large logarithms, as is 
always the case for processes with two or more disparate energy scales. Such logarithms are of 
practical relevance, since they can spoil the convergence of perturbation theory, and they may need 
to be resummed; on the other hand, they are also theoretically interesting, since they are tied to the 
infrared, semiclassical regime: this means that they can be computed in principle to all orders in 
perturbation theory, providing a useful handle on the non-perturbative structure of the theory. 

Studies of the high-energy limit began with non-relativistic potential models, and the first 
powerful tool that was brought to bear was the analytic continuation of amplitudes in the complex 
angular momentum plane [4]. The starting point for these studies is the well-known expansion of 
scattering amplitudes in partial waves. In the case of a four-point amplitude, one writes 

oo 

A{s,t) = 167r£(2/ + l)a/(5)P/(cos0) , (1.1) 
/=o 

where Pi are Legendre polynomials, d is the center-of-mass scattering angle, in terms of which 
t = —s{\ — cos0)/2, and ai{s) are the ^-channel partial wave amplitudes. Using the causality and 
analyticity of the S matrix as fundamental postulates, one can analytically continue the amplitude 
from the physical region {s > and f < 0) to the crossed f-channel, and construct integral repre- 
sentations of the crossed-channel partial-wave amplitudes ai{t). One finds then that the analytic 
properties of di{t), as functions of the angular momentum variable /, have powerful direct impli- 
cations on the high-energy behavior of the original amplitude A{s,t). For example, assuming that 
the only singularities of di{t) in the angular momentum plane are isolated poles, one finds that the 
large-5 asymptotic behavior of the amplitude is given by a simple power law, 

^'^''^^T^) A{s,t) ^ ms'^^'l (1.2) 

where 05(f) is the Regge trajectory. 

This general result follows just from physically motivated postulates about the S matrix, with- 
out reference to an underlying lagrangian field theory, and without resorting to a perturbative expan- 
sion. When one considers a specific field theory model, and uses the tools of perturbation theory, 
one may verify that the structure schematically described by eq. (1.2) does indeed emerge. In a 
gauge theory, the high-energy limit is dominated at tree level by the f-channel exchange of mass- 
less gauge bosons, and it is possible to show, at least at Leading Logarithmic (LL) accuracy [5], 
that virtual corrections dress the f-channel gauge boson propagator according to 




(1.3) 



where the Regge trajectory a{t) is now expressed as a perturbative expansion. A typical example 
of this structure is the resummed expression for the four-gluon amplitude in QCD, which can be 



2 



Infrared singularities at high-energy 



Lorenzo Magnea 



written as 



Mss-^ss A 



28^ 



a(t) 



(1.4) 



ki,Y. The factorized form of eq. (1.4), 



where we label momenta so that s = (^1+^2) and t = {k\— ^3 
which has been shown to hold at NLL accuracy [6] for the real part of the amplitude, is consistent 
with eq. (1.2), and thus with the assumptions that the only singularities in the complex angular 
momentum plane should be isolated poles. The functions Cxjij{ki,kj) are the gluon impact factors, 
and are universal at least to the stated logarithmic accuracy: it is expected that one may use the 
same impact factor in other high-energy processes involving gluons as well as other energetic 
particles, for example quarks. In perturbative QCD, of course, virtual corrections to the amplitude 
ai^e infrared divergent, and these divergences will appear- in both the impact factors and the Regge 
trajectory. To organize these divergences, it is useful employ dimensional regularization (with 
d = A — 2e and e < 0), and express the trajectory as a perturbative series in powers of the d- 
dimensional running coupling. One writes then 



a{t) 



4k 



4n 



a 



(2) + ^(q3) 



(1.5) 



where the J-dimensional running coupling [7] is given by 



a,{-t,e) 



2\ e 



a.(M')+^(a2). 



(1.6) 



The first two perturbative coefficients a^'^ and a^^-* are known [8, 9, 10, 11, 12], and can be written 



as 



wi) 



aP) = Ca 



where bo = (UCa — 2nj)/3, and y^^ are the perturbative coefficients of the cusp anomalous di- 
mension [13], with the overall Casimir eigenvalue Ca scaled out. 

One must naturally wonder to what extent the factorization in eq. (1.4) can be trusted to higher 
logarithmic accuracy, a question which is intimately connected to the underlying assumption that 
the only singularities in angular momentum should be isolated poles. This question was studied in 
the early days of Regge theory, both on general grounds and in specific models (see, for example, 
[2]). In general, one may expect to find cuts as well as poles in the complex angular momentum 
plane. It can be shown that, in a perturbative expansion, these 'Regge' cuts require the f-channel 
exchange of at least two particles, and further require a non-planai^ diagrammatic structure. The 
simplest class of diagrams with all the required features, which appears starting at thr ee loops, is the 
'Mandelstam double-cross', portrayed in Fig. 1. As we shall see, our approach to the high-energy 
limit will lead us to prove that the simple factorized form of the amplitude given in eq. (1.4) must 
break down at NNLL, and the first energy logarithm not correctly predicted by a naive Regge- 
pole-based factorization arises at three loops, from non-planar diagrams, in agreement with the 
expectation that Regge cuts should play a role at precisely that level of accuracy. 



(1.7) 
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Figure 1: The Mandelstam double-cross diagram. 

In the following, we will briefly describe our knowledge of the general structure of infrared 
poles in multi-particle fixed-angle massless gauge theory amplitudes, and we will discuss how this 
knowledge can be applied to the high-energy limit, summarizing the results of [14]. We will show 
that infrared divergences naturally 'Reggeize', for general f -channel color exchanges, and we will 
prove that the simplest form of Reggeization breaks down at NNLL accuracy. Interestingly, the 
general results of Regge theory also have implications on the structure of infrared poles, provid- 
ing non-trivial constraints on possible corrections to the simplest ansatz for the soft anomalous 
dimension, the dipole formula [15, 16]. 

2. The infrared structure of gauge amplitudes 

Let us consider a scattering amplitude in a massless gauge theory, involving n colored particles 
in arbitrary representations of the gauge group. Such an amplitude can be seen as a vector in the 
vector space spanned by all possible color tensors connecting the various representations involved. 
One writes 

Ma,...a„ (^,a.(M'),£) = £M,(^^,a,(At2),£) (q)„....„„, (2.1) 

where the color tensors c j form a basis in the relevant vector space, p, are the particle momenta, 
and a, their color indices. 

Such on-shell matrix elements are plagued by soft and coUinear divergences order by order in 
perturbation theory. Studies performed over the space of several decades (see [17] for a review), 
and summarized in this context in Refs. [15, 18], have shown that infrared singularities factorize. 
Denoting by M the vector with components Mj, one can write [15, 16] the matrix equation 

M(^^,a,(At2),£) =z(^^,a,(At)),£) H [^,^M^''),e^ , (2.2) 

where H is a vector of matching coefficients, which are finite as £ — )• 0, while all infrared singu- 
larities are collected in the matrix factor Z, and /X/ is a factorization scale. A key property of the 
infrared factor Z is that it obeys a (matrix) renormalization group equation with a finite anomalous 
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dimension T, which reads 



lij^zl^M^^),e^ = -z(^,a,(M2),£) r(|,a,(Ai')) . (2.3) 



This equation is easily, if formally, solved in exponential form, as 



Z( ^,a,(Ai^),£ ) =Pexp 



1 f^'dl^^fpi 



A2 \X 



(2.4) 



where we have used as a boundary condition the fact that the tZ-dimensional running coupling 
vanishes in the infrared for d > 4, so that one can set Z(/i = 0) = 1. The anomalous dimension 
matrix F is the centerpiece dictating the infrared structure of multi-particle amplitudes. The key 
result of Refs. [15, 16] is that F itself, in the case of massless particles, obeys a set of exact evolution 
equations which strongly constrain its kinematic dependence, tying it to the color structure of the 
process. The simplest all-order solution to this set of equations is the dipole ansatz for the matrix 
F, which reads 

rdip(^,a.v(A2)) = ^7^(«s.(A2)) £inj^^^x,.T,.-£T^.(a,(A2)) . (2.5) 

Equation (2.5) must be understood as a color operator acting on the hard coefficients H by means 
of the gluon insertion operators T,-, which act exclusively on the color indices of the i-th parti- 
cle, modifying the color structure according to the Feynman rules for single gluon emission [19]. 
The solution is derived under the assumption that the cusp anomalous dimension corresponding 

(R) 

to Wilson lines in color representation R, yi (a^), is proportional, to all orders, to the quadratic 

^ (R) 

Casimir eigenvalue Cr, one may then define the universal function (cCs) according to Yk (0J.v) = 
CrYk (05i)- The functions Yi{as), on the other hand, carry no color structure, and are simply given by 
the anomalous dimensions of the fields corresponding to the the each hard external particle. Equa- 
tion (2.5) has several remarkable features, which cannot be reviewed in detail here. Let us simply 
note that Fjip contains only two-particle coiTclations, which mirrors the structure of one-loop cor- 
rections; it thus preserves the simplicity of single-gluon exchange. This implies vast cancellations 
in higher-order diagrammatic calculations, and means in particular that the color structure Fjip is 
fixed at one loop, so that the path-ordering symbol in eq. (2.4) becomes superfluous. One also sees 
how the structure of double, soft-collinear poles is generated, even as one is simply solving a renor- 
malization group equation, which would ordinarily generate only single poles: indeed, integrating 
the J-dimensional running coupling, accompanied by an extra logarithm of the scale as in eq. (2.5), 
generates precisely such poles. 

Corrections to the dipole ansatz can arise from precisely two sources. First, it is possible that 
the cusp anomalous dimension may cease to be proportional to the quadratic Casimir eigenvalue, 
starting at some loop order. This can in principle happen starting at four loops, where quartic 
Casimirs can arise, however arguments have been given in [16, 21] indicating that this is in fact not 
the case, so that the first correction of this kind should be confined to even higher orders. The only 
other possible corrections to eq. (2.5), which can arise starting at three loops and for at least four 
external particles, must take the form of functions of conformal invariant cross-ratios of external 
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momenta, Piju = {pi ■ pjPk ■ Pi)/{Pi ■ PkPj ■ Pj)- One writes then in general 

;f,«.(A^))^r.,(f 



(2.6) 



The function A has been studied in detail in [20, 21]. A number of constraints can be imposed 
on A, but while they severely limit the available functional forms, they are not sufficient to prove 
that it vanishes, even at the three loop level. Whether such quadrupole corrections do indeed arise, 
starting at three loops, is a question that will probably need to be answered by a direct calculation. 

3. Infrared constraints on high-energy logarithms 

Studying the high-energy limit of the dipole formula is clearly of interest: indeed, the per- 
turbative Regge trajectory is infrared divergent, and one may expect to be able to make general 
predictions, at least for divergent terms, to all orders in perturbation theory. In order to proceed, 
one must first note that the high-energy limit is formally outside the domain of applicability of the 
factorization theorem expressed by eq. (2.2), which was derived for fixed-angle amplitudes, i.e. 
assuming that all kinematic invariants are of the same order, and all large with respect to the strong 
interaction scale A. This objection does not affect our arguments: we can start with a fixed-angle 
configuration, where the factorized expression applies, and continuously deform the kinematics to- 
wards the high-energy limit, while staying in the perturbative regime; what happens is that energy 
logarithms, ln(i^/|f|), become large, and they fail to properly factorize, but no new infrared poles 
are generated. We are thus in a position to make predictions, at least for the infrared divergent part 
of the various ingredients entering the Regge-factorized form of the amplitude. 

Let us consider, for simplicity, the four-point amplitude (analogous considerations apply to 
the multipaiticle amplitudes in multi-Regge kinematics, as shown explicitly in [14]). The key 
property of the infrared factor Z in the Regge limit \t\ is that, to leading power in \t\/s, but to 
all logarithmic accuracies, it factorizes into the product of two factors, one of which carries the s 
dependence and the non-trivial color information, while the second one depends only on t and is 
proportional to the identity in color space. We write this result as 



Z(j,a,(ju2),e) Zi ^-^,a,(At^),e 



(3.1) 



The color-singlet factor Zi contains collinear (as well as soft) divergences to be associated with 
impact factors in a Regge-factorized expression. The central piece of eq. (3.1) is the 5-dependent 
matrix Z, which can be written in a very general and elegant way as 



Z y-,as{}JL ),£) = exp { K[ a.v(/x ),£ 



In 



(3.2) 



Here we have introduced notations for the color insertion operators coiTcsponding to the color 
representations exchanged in the t and s channel [22], defining 



T, = Ti+T2 , T, = Ti+T3, 



(3.3) 
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for the scattering process 1+2 — > 3+4. Furthermore, we have introduced the notation 



K{ 



(3.4) 



It is straightforward to show that eqs. (3.1) and (3.2) imply LL Reggeization of infrared poles in 
?-channel exchanges for the scattering of generic color representations. Indeed, at LL accuracy, 
one may neglect the 5-channel phase in the exponent of eq. (3.2), which starts contributing at NLL. 
One can then write the full amplitude as 



M( —,a,{ix^),e 



^ = exp|^(a,(jU^),e 



In 



(3.5) 



Whenever the Born amplitude (which gives the leading-order contribution to H) is dominated by 
^-channel exchanges at leading power in l^j/i, each such exchange is an eigenstate of the ^-channel 
color operator Tj, with an eigenvalue given by the appropriate Casimir operator. For gluon-gluon 
scattering, for example, the hard scattering is dominated in the high-energy limit by gluon ex- 
change, and one may write 



CaH^ 



gg^gg 



+ 



\t\/s). 



(3.6) 



This automatically leads to Reggeization, yielding 

M^g^Sg 



(3.7) 



and thus providing an all-order expression for the Regge trajectory, given by eq. (3.4), reproducing 
early results [12]. 

Since the factorization in eq. (3.2) is valid at leading power in \t\/s, we can in principle extend 
our predictions to all subleading energy logarithms. One way to do it is to use the Baker-Campbell- 
Haussdorff formula to express Z as a product of exponentials. One finds 



Z( -,a,(At^),£ 



KTr 



exp {iTi/r T^} exp 



-i-K^ In 
2 



X exp 



27iHn( — ] \tI[tIt 



+ iTT In^ 



T?, [T: 



t 



fTi2 



(3.8) 



. Each exponential factor in eq. (3.8) carries increasing powers of the function K{as.,e) ~ (Xs/e. 
One observes that, as expected, only the first factor contributes at LL. At NLL accuracy, one finds a 
series of commutators, which in general are not diagonal in the f-channel basis, and will thus break 
the naive form of Regge factorization expressed for example by eq. (1.4). All these commutators 
however appear as phases, so that they will contribute only to the imaginary part of the amplitude: 
we see, on the other hand, that infrared poles for the real part of the amplitude Reggeize at NLL 
accuracy for ai^bitraiy color exchanges in the t channel, extending the results of [6]. At NNLL 
accuracy, we predict that the simple form of Regge factorization, which is based on the assumption 
that only Regge poles arise in the angular momentum complex plane, must break down also for the 
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real part of the amplitude. The first color operator effecting this breakdown arises at iff{a^), and it 
is given by 

(3.9) 



The most important thing to note concerning eq. (3.9) is that Regge-breaking terms start arising at 
three loops, at NNLL level, and they involve non-planar contributions to the amplitude, as implied 
by the presence of color commutators; furthermore, they arise from the squaring of the ^-channel 
phase in eq. (3.2). All these features closely match what is expected from the presence of Regge 
cuts in the angular momentum plane, as discussed in Sec. 1. A second important point is that 
possible corrections to the dipole formula arising at three loops (the function A in eq. (2.6)) can- 
not rescue naive Regge factorization, since they contribute terms of order /e to the amplitude, 
whereas eq. (3.9) contains terms of order cc^ /s^. We also note that a trace of this Regge-breaking 
mechanism must arise already at two loops in terms with no energy logarithms, which are gener- 
ated for example from the expansion of the second factor in eq. (3.8). Terms of precisely this form, 
proportional to ajn^/e^, where indeed detected in [10] by direct calculation. Finally, we note that 
double color commutators terms such as those generated by eq. (3.9) have recently been shown to 
play a role [24] in the breakdown of strict coUinear factorization described in [25] : in that context, 
these terms arise from the contributions of Glauber gluons which survive in finite hadronic cross 
sections which are not fully inclusive, in the form of super-leading logarithms [26]. 



4. High-energy constraints on infrared poles 

We close by briefly illustrating how the known structure of the Regge limit can be used to im- 
pose further constraints on quadrupole corrections to the dipole formula, which may arise starting 
at the three loop order. Considering again for simplicity the four-point amplitude, and following 
Ref. [20], we note that the thi^ee relevant conformal cross ratio can be written as 



Pl234 = 




(4.1) 



where the appropriate phases have been kept, and the constraint P1234P1342P1423 = 1 is verified. 
Ref. [20], after examining the constraints imposed on the function A in eq. (2.6) by soft-coUineai^ 
factorization, Bose symmetry and transcendentality, found a small set of functions that would still 
give allowed contributions. The functions considered in [20] have simple dependences on loga- 
rithms or polylogarithms of the conformal cross ratios in eq. (4.1), and one can study their high- 
energy limit. Defining L = \n{s/\t\), and Ljji^i = In pijj^i, one readily sees that, to leading power in 
s/\t\, the logarithms Ljjti become 

L1234 = 2(L-i7r), L1342 = -2L, Li423 = 2i7r. (4.2) 

One can now examine what happens to the non-trivial examples of A that were discussed in [20]. 
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The simplest case, involving only logarithms, is the function 



jadejche ^2^^^ / ^^^^^ ^2^^^ _|_ ^2^^^ ^^^^^ | ^4 j-j 



+ J J i-1423 (^^1234^-1342 + ^1234^1342 
+ f^'^T'^'' 1-\m2 (^1423^1234 +^1423^1234^ 

It is now immediate to verify that, in the high-energy limit 

+ (linL^ - 3k^L^ - in^L) /™7'"«^j +ff(\t/s\). (4.4) 
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The known structure of the Regge limit forbids the appearance of such a function, since it would 
generate super-leading high-energy logarithms of the form a"L"+^ at three loops and beyond. 
Studying the high-energy limit is sufficient to rule out, individually, all explicit examples of A 
that were discussed in [20]. It is clear however that it may still be possible to consider linear com- 
binations of the same functions, constructed precisely in order to cancel the logarithms that happen 
to be constrained by Regge theory arguments. Indeed, a recent study [21] has identified sets of 
functions that explicitly satisfy all known constraints, including those arising from the high-energy 
limit. The question of the existence of non-vanishing corrections to the dipole formula at three 
loops and beyond remains therefore open. 
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